Back

Computer Methods in Biomechanics and Biomedical Engineering

Informa UK Limited

All preprints, ranked by how well they match Computer Methods in Biomechanics and Biomedical Engineering's content profile, based on 10 papers previously published here. The average preprint has a 0.01% match score for this journal, so anything above that is already an above-average fit. Older preprints may already have been published elsewhere.

1
Investigating the Role of Neck Muscle Activation and Neck Damping Characteristics in Brain Injury Mechanism

Bahreinizad, H.; Chowdhury, S.

2024-01-26 bioengineering 10.1101/2023.11.15.567289 medRxiv
Top 0.1%
6.8%
Show abstract

PurposeThis study aimed to investigate the role of neck muscle activity and neck damping characteristics in traumatic brain injury (TBI) mechanisms. MethodsWe used a previously validated head-neck finite element (FE) model that incorporates various components such as scalp, skull, cerebrospinal fluid, brain, muscles, ligaments, cervical vertebrae, and intervertebral discs. Impact scenarios included a Golf ball impact, NBDL linear acceleration, and Zhangs linear and rotational accelerations. Three muscle activation strategies (no-activation, low-to-medium, and high activation levels) and two neck damping levels by perturbing intervertebral disc properties (high: hyper-viscoelastic and low: hyper-elastic) strategies were examined. We employed Head Injury Criterion (HIC), Brain Injury Criterion (BrIC), and maximum principal strain (MPS) as TBI measures. ResultsIncreased neck muscle activation consistently reduced the values of all TBI measures in Golf ball impact (HIC: 4%-7%, BrIC: 11%-25%, and MPS (occipital): 27%-50%) and NBDL study (HIC: 64%-69%, BrIC: 3%-9%, and MPS (occipital): 6%-19%) simulations. In Zhangs study, TBI metric values decreased with the increased muscle activation from no-activation to low-to-medium (HIC: 74%-83%, BrIC: 27%-27%, and MPS (occipital): 60%-90%) and then drastically increased with further increases to the high activation level (HIC: 288%-507%, BrIC: 1%-25%, and MPS (occipital): 23%-305%). Neck damping changes from low to high decreased all values of TBI metrics, particularly in Zhangs study (up to 40% reductions). ConclusionOur results underscore the pivotal role of neck muscle activation and neck damping in TBI mitigation and holds promise to advance effective TBI prevention and protection strategies for diverse applications.

2
Strange effects of activation dynamics on musculoskeletal trajectory optimization

van den Bogert, A. J.

2025-02-05 bioengineering 10.1101/2025.01.30.635759 medRxiv
Top 0.1%
6.8%
Show abstract

Muscle activation dynamics is usually described with a nonlinear differential equation, such that the activation occurs faster than deactivation. When a muscle excitation input switches rapidly between two values, this produces a "pumping" effect that raises the activation above the mean excitation input. This effect will favor rapid switching strategies when excitation is used in the optimization objective for human movement. The optimal strategy has infinitely fast switching, and is difficult to obtain with direct collocation methods. This problem can be eliminated by using activation, rather than excitation, in the optimization objective.

3
Astrocytes in white matter respond to tensile cues during cortical folding: a numerical study

Taneja, K.; Saito, K.; Kawasaki, H.; Holland, M. A.

2025-10-19 bioengineering 10.1101/2025.10.17.683172 medRxiv
Top 0.1%
6.7%
Show abstract

Our understanding of the process of formation of gyri (ridges) and sulci (furrows) in the cerebrum during development is moving beyond the role of neurons. Glial cells such as astrocytes, which are the most common cell type in the brain, are especially prominent under the gyri and have been shown to be essential for gyrification in ferrets. Their dysfunction has been linked to a host of neurodevelopmental diseases and disorders in humans, leading to abnormal folding patterns, hence, it is crucial to understand their role in the mechanics of folding. In this work, we propose two hypotheses of how astrocytes affect cortical folding. Our previous study demonstrated that astrocytes proliferate in the white matter (subcortex), and that inhibiting this process impairs gyrification. This leads to the pushing hypothesis, where astrocytes push up the cortex outwards, leading to formation of folds. On the other hand, ex vivo studies demonstrate areas of the cortex and subcortex that experience tension, due to the differential growth between materials in the brain tissue. This leads to the pulling hypothesis, where astrocytes experience tension from the surrounding tissue, leading to their proliferation, distribution of tensile stresses, and initiation of growth in those regions. Using the theory of finite growth, we implement these hypotheses via morphogenetic growth (pushing) and stress-driven (pulling) growth criteria. We find that morphological trends during development between the pushing and pulling effects are not dissimilar, with a comparable gyrification index. The stress distributions from both models also show common features, but the pulling effect shows tension in the subcortex, which matches trends observed in experiments. Therefore, it is more likely that the astrocytes affect gyrification by proliferating as a response to tensile cues and decrease the tension they experience, leading to deeper folds, rather than astrocytes independently proliferating under a gyri and then pushing the cortex up.

4
COVID-19 isolation and containment strategies for ships: Lessons from the USS Theodore Roosevelt outbreak

Stoddard, M.; Johnson, K.; White, D.; Nolan, R.; Hochberg, N.; Chakravarty, A.

2020-11-07 epidemiology 10.1101/2020.11.05.20226712 medRxiv
Top 0.1%
6.7%
Show abstract

The control of shipborne disease outbreaks represents a vexing but entirely predictable challenge at the start of any pandemic. Passenger ships, with large numbers of people confined in close quarters, can serve as incubators of disease, seeding the pandemic across the globe as infected passengers return home. Short-term steps taken by local authorities can exacerbate this problem, creating humanitarian crises and worsening the scale of the outbreak. In this work, we have undertaken a model-based examination of the USS Theodore Roosevelt outbreak to understand the dynamics of COVID-19 spread aboard the aircraft carrier. We have used a series of counterfactual "what-if" analyses to better understand the options available to public health authorities in such situations. The models suggest that rapid mass evacuation and widespread surveillance testing can be effective in these settings. Our results lead to a set of generalizable recommendations for disease control that are broadly applicable to the current COVID-19 crisis as well as to future pandemics.

5
The effects of muscle fibre type distribution on gait biomechanics: A predictive simulation study

Daehlin, T. E.; Ross, S. A.; De Groote, F.; Wakeling, J. M.

2026-04-15 bioengineering 10.64898/2026.04.13.718234 medRxiv
Top 0.1%
6.1%
Show abstract

AO_SCPLOWBSTRACTC_SCPLOWMuscle fibre type distribution influences both the metabolic and contractile properties of individual muscles. However, as humans tend to self-optimize their gait pattern to minimize cost of transport, these changes in muscle properties may influence gait biomechanics in manners that are difficult to isolate in in vivo experiments. The purpose of this study was to predict the influence of muscle fibre type distribution on the metabolic cost and biomechanics of simulated walking and running. We implemented a muscle model that could predict recruitment of slow and fast twitch muscle fibres in a framework for predictive musculoskeletal simulation. Subsequently, we employed the framework to investigate how metabolic cost of transport, stride length, stride frequency, and mechanical work performed by slow and fast twich muscle fibres were influenced by fibre type distribution across locomotion speeds from 1.0 to 4.5 m {middle dot} s-1. Our results predict that cost of transport increases as slow twitch area fraction decreases, while stride length and frequency was minimally affected by fibre type distribution at speeds resulting in walking. In contrast, fibre type distribution interacts with locomotion speed at speeds resulting in running. Specifically, we predict the existence of a threshold speed below which cost of transport decreases with an increasing proportion of slow twitch fibres, while cost of transport increases with increasing proportions of slow twitch fibres above it. The shift in fibre type distribution was accompanied by an increase in stride frequency and decrease in stride length. These shifts in spatiotemporal characteristics appear to allow the muscles to operate at speeds close to those that achieve peak mechanical efficiency. Taken together, the results of this study predict that muscle fibre type distribution may influence both the energetics and biomechanics of gait, and that this influence is dependent upon the locomotion speed.

6
Vehicle safety tests, rankings, curb weight, and fatal crash rates: automatic emergency brakes associated with increased death rates

Robertson, L. S.

2022-12-12 epidemiology 10.1101/2022.12.08.22283253 medRxiv
Top 0.1%
5.7%
Show abstract

The European New Car Assessment Programme (Euro NCAP), the Insurance Institute for Highway Safety (IIHS), and the U.S. National Highway Traffic Safety Administration (NHTSA) each publish safety ratings of new passenger vehicles based on crash test results and crash avoidance technology soon after they are introduced into the market. IIHS alone singles out vehicles for "Top Safety Pick". The Institute also periodically lists driver death rates of vehicles during their first few years of use, accompanied by assertions that larger, heavier passenger vehicles are safer. Median death rates by vehicle size in the Institutes data do not support the assertion. This study examines the association of vehicle weight to the risk of all deaths in fatal crashes where specific makes and models were involved, controlling statistically for the Institutes vehicle safety ratings as well as NHTSA ratings of full-frontal crash tests and rollover propensity, lane retention warnings, adaptive cruise control, and automatic braking. Increased weight is slightly related to a lower fatal risk in the 2014-2017 models but not the 2018-2019 models. Lane retention warnings and IIHSs higher ratings of crashworthiness are related to reduced death risk. Adaptive cruise control is not associated with fatal crash risk. Automatic emergency braking (AEB) technology rated "superior" is a criterion for an IIHS "Top Safety Pick" but it is correlated to higher fatal crash involvement of vehicles that have the technology as optional or standard equipment. IIHS tests of AEB systems are conducted only at low speeds. The statistical results are not definitive but suggest that a rating system based on modeling death risk prediction from more detailed data from various tests such as Euro NCAP tests of AEB systems could better inform consumer choice of vehicles.

7
Maize Stem Buckling Failure is Dominated by Morphological Factors

Stubbs, C. J.; Larson, R.; Cook, D. D.

2019-11-07 bioengineering 10.1101/833863 medRxiv
Top 0.1%
5.6%
Show abstract

The maize (Zea mays) stem is a biological structure that must both balance biotic and structural load bearing duties. These competing requirements are particularly relevant in the design of new bioenergy crops. With the right balance between structural and biological activities, it may be possible to design crops that are high-yielding and have digestible biomass. But increased stem digestibility is typically associated with a lower structural strength and higher propensity for lodging. This study investigates the hypothesis that geometric factors are much more influential in determining structural strength than tissue properties. To study these influences, both physical and in silico experiments were used. First, maize stems were tested in three-point bending. Specimen-specific finite element models were created based on x-ray computed tomography scans. Models were validated by comparison with in vitro data. As hypothesized, geometry was found to have a much stronger influence on structural stability than material properties. This information reinforces the notion that deficiencies in tissue strength could be offset by manipulation of stalk morphology, thus allowing the creation of stalks with are both resilient and digestible. HighlightThis study utilized physical and in silico experiments to confirm that geometric parameters are far more influential in determining stalk strength than mechanical tissue stiffnesses.

8
Prediction of Covid-19 Infections Through December 2020 for 10 US States Incorporating Outdoor Temperature and School Re-Opening Effects-August Update

Newell, T. A.

2020-09-18 epidemiology 10.1101/2020.09.14.20193821 medRxiv
Top 0.1%
5.5%
Show abstract

End-of-August updates for Covid-19 infection case predictions for 10 US States (NY, WA, GA, IL, MN, FL, OH, MI, CA, and NC) are compared to actual data. Several states that experienced significant summer surges gained control of accelerating infection spread during August. The US as a whole and the 10 States investigated continue to follow periods of linear infection growth that defines a boundary separating accelerated infection growth and infection decay. August 31 predictions (initiated July 27, 2020) for four States (NY, WA, MI and MN) are within 10% of actual data. Predictions for four other States (GA, IL, CA, and OH) agree between 10 and 20% of actual data. Predictions for two States (FL and NC) are greater than 20% different from actual data. Systematic differences between predictions and actual data are related to the impact of the June-July summer surge, and human behavior reactions (ie, increased face mask usage and distancing) to accelerated infection growth. Outdoor temperature effects and school re-opening effects are not apparent nor expected for August. Human behavior parameters (Social Distance Index, SDI, and disease transmission efficiency, G, parameters) are adjusted to mirror August data. Comparisons of actual versus predicted daily new infection cases display the complexity of SARS-CoV-2 transmission.

9
The role of falx cerebri in the selective vulnerability of splenium within the corpus callosum

zhou, z.; kleiven, s.

2026-06-02 bioengineering 10.64898/2026.05.31.729036 medRxiv
Top 0.1%
5.5%
Show abstract

The corpus callosum is the largest white matter structure connecting the two cerebral hemispheres and is anatomically divided into three major subregions along the anteroposterior axis: the genu, midbody, and splenium. The splenium is frequently affected in traumatic head impacts, yet the biomechanical basis for this selective vulnerability remains poorly understood. Clinical studies have long hypothesized that the falx cerebri contributes to the splenial susceptibility because of its close anatomical relationship with the posterior corpus callosum, although direct verification is lacking. To address this, a high-resolution finite element head model with explicit representations of the genu, midbody, and splenium was employed. Two model variants, differing only in the presence or absence of an anatomically and mechanically detailed falx, were used to simulate ten head impacts covering a range of loading directions and severities. Peak strain, strain rate, and shear stress were quantified in each corpus callosum subregion and compared using linear mixed-effects models. The results showed that inclusion of the falx altered the regional distribution of mechanical responses within the corpus callosum. Across the simulated impacts, the splenium consistently exhibited greater strain, strain rate, and shear stress than the genu and midbody when the falx was present. In contrast, these preferentially larger splenial deformation were not consistently observed when the falx was absent. Statistical analyses demonstrated significant region-dependent effects of the falx, with falx-induced increases in strain, strain rate, and shear stress being significantly greater in the splenium than in the genu and midbody (p < 0.05). These findings verified the hypothesis that the falx selectively amplified mechanical loading within the splenium, thereby contributing to its heightened vulnerability to injury. This work provides a plausible biomechanical explanation for the frequent involvement of the splenium in brain trauma patients and highlights the heterogeneous influence of the falx on mechanical responses across corpus callosum subregions.

10
Estimating the scale of COVID-19 Epidemic in the United States: Simulations Based on Air Traffic directly from Wuhan, China

Li, D.; Lv, J.; Botwin, G.; Braun, J.; Cao, W.; Li, L.; McGovern, D. P. B.

2020-03-08 epidemiology 10.1101/2020.03.06.20031880 medRxiv
Top 0.1%
5.4%
Show abstract

IntroductionCoronavirus Disease 2019 (COVID-19) infection has been characterized by rapid spread and unusually large case clusters. It is important to have an estimate of the current state of COVID-19 epidemic in the U.S. to help develop informed public health strategies. MethodsWe estimated the potential scale of the COVID-19 epidemic (as of 03/01/2020) in the U.S. from cases imported directly from Wuhan area. We used simulations based on transmission dynamics parameters estimated from previous studies and air traffic data from Wuhan to the U.S and deliberately built our model based on conservative assumptions. Detection and quarantine of individual COVID-19 cases in the U.S before 03/01/2020 were also taken into account. A SEIR model was used to simulate the growth of the number of infected individuals in Wuhan area and in the U.S. ResultsWith the most likely model, we estimated that there would be 9,484 infected cases (90%CI 2,054-24,241) as of 03/01/2020 if no successful intervention procedure had been taken to reduce the transmissibility in unidentified cases. Assuming current preventive procedures have reduced 25% of the transmissibility in unidentified cases, the number of infected cases would be 1,043 (90%CI 107-2,474). ConclusionOur research indicates that, as of 03/01/2020., it is likely that there are already thousands of individuals in the US infected with SARS-CoV-2. Our model is dynamic and is available to the research community to further evaluate as the situation becomes clearer.

11
Modelling Human Gait using a Nonlinear Differential Equation

Schmalz, J.; Paul, D.; Shorter, K.; Cooper, M.; Murphy, A.

2021-03-17 bioengineering 10.1101/2021.03.16.435713 medRxiv
Top 0.1%
5.4%
Show abstract

We introduce an innovative method for the investigation of human gait, which is based on the visualisation of the vertical component of the movement of the centre of mass during walking or running, in the space of the coordinates position, velocity, and acceleration of the centre of mass. Collected data has been numerically approximated by the best fitting curve for a non-linear model. The resulting equation for the best fitting plane or curve in this space is a differential equation of second order. The model that we suggest is a Duffing equation with coefficients that depend on the height of a walker or runner and on the angular frequency of the oscillation. Statistics about the distribution of the Duffing stiffness depending on the speed is presented. 1 Author SummaryWe study the human gait modelled by the movement of the centre of mass of the test person. This is an example of a biological process which can be considered as a periodical dynamic system. Roughly, this movement behaves in a similar way to a vibrating mass suspended on a spring, but it is more complex. The vertical component of the movement during walking or running can be visualised as an oscillogram: a graph of the position as a function of time. We suggest a visualisation of the data in 3D space, where the coordinates describe position, velocity, and acceleration. Our new visualisation method allows us to model the movement of a persons centre of mass by a nonlinear differential equation. The resulting curve for an ideal spring-mass movement, without viscosity or external force, is an ellipse in the suggested 3D space. The shape of the data curve shows at which position an additional force was applied, or the movement slowed down. Some deviations are common for all test persons and others are different. In the future we plan to investigate the reasons for these deviations, such as different running techniques or the presence of injuries.

12
The muscle coordination required for efficient locomotion scales with body size

Latreche, A.; Ross, S. A.; Dick, T. J. M.; Konow, N.; Biewener, A. A.; Wakeling, J. M.

2026-05-03 bioengineering 10.64898/2026.04.30.722018 medRxiv
Top 0.1%
5.3%
Show abstract

AO_SCPLOWBSTRACTC_SCPLOWMuscle efficiency decreases with increasing size, largely due to a relative decrease in its mechanical output. Muscle mechanical output depends on its activation, strain, and strain rate and thus varies between different muscles within a limb during locomotion. Distinct muscle coordination patterns are required for efficient cycling, and so we would expect that the coordination patterns for efficient cycling or indeed locomotion would change across animal sizes. We tested whether muscle coordination would change with muscle size using data derived from human cycling: this paradigm allowed for controlled changes in both crank torque and cadence, allowing the multifactorial problem of muscle power output to be decomposed. We used kinematic and pedal data from 12 cyclists undergoing steady pedalling at cadences from 80 to 140 r.p.m. and generated musculoskeletal simulations of their movements. We introduced novel multisegment muscle models in the simulation that incorporated the internal muscle mass and thus accounted for the scaling effects of muscle tissue inertia. We solved the simulations for the muscle activity that was required to minimise the metabolic cost during cycling for each condition. The masses of the muscle models were scaled across five orders of magnitude. The predicted muscle activations were classified by Principal Component analysis to identify whether the coordination of muscle activity was modulated across models with different sized muscles. Analysis of variance revealed significant changes in coordination at the large-scale factors. This study shows how the coordination of muscle activity during locomotion will likely change across a range of body sizes due to the non-linear effects of the inertial mass within the muscle tissues.

13
A Numerical Method to Compute Brain Injury Associated with Concussion

Bastien, C.; Scattina, A.; Neal-Sturgess, C. E.; Panno, R.; Shrinivas, V.

2022-10-27 bioengineering 10.1101/2022.10.26.513868 medRxiv
Top 0.1%
5.2%
Show abstract

Concussion can result from various events in everyday life, including falls, sports collisions, and motor vehicle accidents, which could lead to the disruption of neuronal cell membranes and axonal stretching, leading to a neuro-metabolic cascade of molecular changes in the brain. There is currently no agreement on which computational method can assess such low-level injuries. This paper demonstrates for the first time that the Peak Virtual Power (PVP), based on the Clausius-Duhem inequality, assuming that the injury is represented by the irreversible work in a human body, could be a candidate to capture brain distortion related to concussion. The work is based on the evaluation of the PVP via reconstruction of three NFL helmet-to-helmet impacts by means of finite element analysis, using validated Biocore helmet models fitted with calibrated Hybrid III headforms against linear and angular acceleration impact corridors, which were defined as realistic impact conditions for each collision scenario. Once the exact impact parameters were defined, the Hybrid III headform was replaced with a validated THUMS 4.02 human head model in which the PVP was computed for each head at the corpus callosum and midbrain locations. The results indicate that mild and severe concussions could be prevented for lateral collisions and frontal impacts with PVP values lower than 0.928mW and 9.405mW, respectively, and no concussion would happen in the head vertical impact direction for a PVP value of less than 1.184mW. This innovative method proposes a new paradigm to improve helmet designs, assess sports injuries and improve peoples wellbeing. HighlightsO_LIPeak Virtual Power method can capture brain distortion related to concussion C_LIO_LIConcussion is extracted from corpus callosum and midbrain locations of THUMS4.02 C_LIO_LIPeak power in midbrain less than 1.184mW for a vertical impact leads to no concussion C_LIO_LIPeak power in midbrain more than 0.928mW for a lateral impact leads to concussion C_LIO_LIPeak power in midbrain more than 9.405mW for a front impact leads to concussion C_LI

14
Using physiologically-based models to predict in vivo skeletal muscle energetics

Konno, R. N.; Lichtwark, G. A.; Dick, T. J.

2024-05-23 bioengineering 10.1101/2024.05.21.595083 medRxiv
Top 0.1%
5.1%
Show abstract

Understanding how muscles use energy is essential for elucidating the role of skeletal muscle in animal locomotion. Yet, experimental measures of in vivo muscle energetics are challenging to obtain, so physiologically-based muscle models are often used to estimate energy use. These predictions of individual muscle energy expenditure are not often compared to indirect whole body measures of energetic cost. Here, we examined and illustrated the capability of physiologically-based muscle models to predict in vivo measures of energy use. To improve model predictions and ensure a physiological basis for model parameters, we refined our model to include data from isolated muscle experiments. Simulations were performed to capture three different experimental protocols, which involved varying contraction frequency, duty cycle, and muscle fascicle length. Our results demonstrated that these models are able capture the general features of whole body energetics across contractile conditions, but tended to under predict the magnitude of energetic cost. Our analysis revealed that when predicting in vivo energetic rates across contractile conditions, the model was most sensitive to the force-velocity parameters and the data informing the energetic rates when predicting in vivo energetic rates across contractile conditions. This work highlights it is the mechanics of skeletal muscle contraction that govern muscle energy use.

15
The Powered Simplest Walking Model Explains the Different Vertical Ground Reaction Force Amplitudes at Elevated Walking Speeds

Hosseini-Yazdi, S.-S.

2024-08-01 bioengineering 10.1101/2024.07.16.603707 medRxiv
Top 0.1%
5.1%
Show abstract

Understanding the vertical ground reaction force (vGRF) profile offers important insight into how humans regulate mechanical work during walking. Although the characteristic double-hump vGRF pattern is well documented, the mechanical factors underlying asymmetry in peak amplitudes and midstance trough timing remain unclear. Using a simple powered walking model and an inverted pendulum simulation with constant hip torque, we examined how step-transition work--collision and push-off--shapes the vGRF trajectory. We further compared these predictions to empirical data spanning walking speeds from 0.8-1.4 m. s-1. The simple walking model predicted symmetric vGRF profiles across speeds because collision and push-off impulses were equal, resulting in passive single-support motion. In contrast, adding hip torque within the pendular model produced stance-phase asymmetries, shifting the vGRF trough earlier when torque added energy and later when torque dissipated energy. Empirical analysis revealed that collision and push-off impulses were generally unequal except at one speed, producing asymmetric vGRF peaks. At low speeds, push-off exceeded collision; at high speeds, the reverse occurred, consistent with a need for compensatory single-support positive work. These mechanical imbalances predicted systematic shifts in trough timing toward the dominant impulse. Therefore, we propose the Vertical GRF Trough Timing Index (vGRF-TTI), combined with collision and push-off peak amplitudes, as a clinically meaningful outcome capturing the balance of step-transition work. Earlier troughs with elevated collision peaks indicate impaired push-off or constrained gait conditions, whereas later troughs with larger push-off peaks reflect compensatory or enhanced propulsion. These metrics provide sensitive, mechanism-based indicators of gait efficiency and neuromotor control.

16
Scaling contact force parameters across body size, limb count, and number of contact spheres

van Bijlert, P. A.

2025-11-29 biophysics 10.1101/2025.11.26.690874 medRxiv
Top 0.1%
5.0%
Show abstract

A popular way to model contact interactions in musculoskeletal simulations uses Hertz theory applied to contact spheres, with Hunt Crossley based dissipation. Suitable contact parameters for dynamic simulations will be highly dependent on the morphology, scale, materials, and movement in question. Inappropriate parameter choices can manifest in unpredictable ways during simulations, potentially resulting in misinterpretations or failed simulations. Here, I demonstrate that both the plane strain modulus and the dissipation parameters are not scale invariant. I derive equations to scale the contact parameters in dimensionless form, which allows accounting for differences in body size, number of legs, contact sphere radius, and number of spheres per foot. As a demonstration of this scaling approach, I scale the contact parameters of a 62 kg human to a 500 kg human, a mouse (0.02 kg), an emu (37.8 kg), a horse (545 kg), and a giraffe (1190 kg), and demonstrate that geometrically and dynamically similar contact behaviour is achieved in all cases. The scaling approach presented here can be used to scale parameters known to work for one model to a completely different model, which is particularly useful in studies that simulate the effects of allometric scaling. I also provide equations to estimate suitable contact parameters for a model directly, without using a different model as a starting point. The limitations of Hertz Hunt Crossley contact models in biomechanical simulations are discussed. Lastly, I derive dimensionless expressions and scaling guidelines for the smoothed contact force implementation "SmoothSphereHalfSpaceForce" in the popular biomechanical simulator OpenSim.

17
Triphasic Thrombosis Model: A Computational Study of Type B Aortic Dissection

Gupta, I.; Schanz, M.; Ricken, T.

2024-07-04 biophysics 10.1101/2024.05.07.592918 medRxiv
Top 0.1%
5.0%
Show abstract

Thrombosis refers to the formation of a thrombus, or a blood clot, within the body, which can occur either partially or completely. It serves as a crucial indicator of the severity of a patients medical condition, with the location and characteristics of thrombosis dictating its clinical implications. Hence, accurate diagnosis and effective management of thrombosis are paramount. In our current investigation, we incorporate the porous attributes of a thrombus using the Theory of Porous Media. This involves dividing the aggregate into solid, liquid, and nutrient phases and utilising volume fractions to capture microstructural details. Fluid flow through the porous media is modelled using a modified Darcy-Brinkman type equation, with interaction terms within balance equations facilitating the modelling of the mass exchange and other phase interactions. The shorter time scales are neglected. We present a comprehensive framework of equations and assumptions governing the behaviour of a strongly coupled multiphasic porous medium problem. Additionally, we introduce scenarios involving type B Aortic Dissection and false lumen geometries, providing a detailed outline of the problem setup. Thereafter, we present the potential of the model for thrombi growth. The simulation results are compared with velocity plots aligning with Magnetic Resonance Imaging data for three distinct cases with varying entry and exit tear sizes. Consequently, our proposed model offers a promising and reasonable approach for numerically simulating thrombosis and gaining insights into the underlying growth mechanics.

18
Did long COVID increase road deaths in the U.S.?

Robertson, L. S.

2023-10-16 epidemiology 10.1101/2023.10.11.23296868 medRxiv
Top 0.1%
5.0%
Show abstract

ObjectiveTo examine data on COVID-19 disease associated with a 10 percent increase in U.S. road deaths from 2020 to 2021 that raises the question of the potential effect of pandemic stress and neurological damage from COVID-19 disease. MethodsPoisson regression was used to estimate the association of recent COVID-19 cases, accumulated cases, maximum temperatures, truck registrations, and gasoline prices with road deaths monthly among U.S. states in 2021. Using the regression coefficients, changes in each risk factor from 2020 to 2021 were used to calculate expected deaths in 2021 if each factor had remained the same as in 2020. ResultsCorrected for the other risk factors, road deaths were associated with accumulated COVID-19 cases but not cases in the previous month. More than 20,700 road deaths were associated with the changes in accumulated COVID-19 cases but were substantially offset by about 19,100 less-than-expected deaths associated with increased gasoline prices. ConclusionsWhile more research is needed, the data are sufficient to warn people with "long COVID" to minimize road use. What is already known about this topicPrevious short-term fluctuations in road deaths are related to changes in temperature, fuel prices, and truck registrations. What this study addsCorrected for other risk factors, the monthly changes in road deaths from 2020 to 2021 in U.S. states were associated with cumulative COVID-19 cases. How this study might affect research, practice, or policyStudies are needed to distinguish the potential relative effects of neurological damage as well as the stress of coping with the pandemic on driving, walking, and bicyclist behavior. Warning people with "long covid" about road risk is warranted.

19
Predicted effects of summer holidays and seasonality on the SARS-Cov-2 epidemic in France

Duchemin, L.; Veber, P.; Paris, M.; Boussau, B.

2020-07-07 epidemiology 10.1101/2020.07.06.20147660 medRxiv
Top 0.1%
4.9%
Show abstract

1The SARS-CoV-2 epidemic in France has had a large death toll. It has not affected all regions similarly, since the death rate can vary several folds between regions where the epidemic has remained at a low level and regions where it got an early burst. The epidemic has been slowed down by a lockdown that lasted for almost eight weeks, and individuals can now move between metropolitan French regions without restriction. In this report we investigate the effect on the epidemic of summer holidays, during which millions of individuals will move between French regions. Additionally, we evaluate the effect of strong or weak seasonality and of several values for the reproduction number on the epidemic, in particular on the timing, the height and the spread of a second wave. To do so, we extend a SEIR model to simulate the effect of summer migrations between regions on the number and distribution of new infections. We find that the model predicts little effect of summer migrations on the epidemic, because the number of migrating infectious individuals are low as a consequence of the lockdown. However, all the reproduction numbers above 1.0 and the seasonality parameters we tried result in a second epidemic wave, with a peak date that can vary between October 2020 and April 2021. If the sanitary measures currently in place manage to keep the reproduction number below 1.0, the second wave will be avoided. If they keep the reproduction number at a low value, for instance at 1.1 as in one of our simulations, the second wave is flattened and could be similar to the first wave.

20
Tail wags the dog is unsupported by biomechanical Modeling of Canidae Tails Use during Terrestrial Motion

Rottier, T.; Schulz, A.; Sohnel, K.; McCarthy, K.; Fischer, M.; Jusufi, A.

2022-12-31 biophysics 10.1101/2022.12.30.522334 medRxiv
Top 0.1%
4.8%
Show abstract

Dogs and other members of Canidae utilize their tails for different purposes, including agile movements, such as running and jumping. In this study, we utilized motion capture biomechanical data of a border collie executing an agile rotational jump maneuver. This data created a 17-segment biomechanical model of the border collies (Canis familiaris) limb movement during agile jumps. This model was verified by comparing it to the biomechanical movement and fitting the dogs agile task with an RMSE less than 2.5%. Using this joint model, we held specific segments constant to view their inertial impact on the dog during the aerial phase of jumping. Results suggest that the tail, hind limbs, and fore limb provides little to no inertial advantage during these rotational jump maneuvers. The tail of dogs likely does have a minimal impact on inertia, the opposite of animals like the gecko. This work could alleviate unknown biomechanical use of the tails to understand the behavioral biomechanics of lesser-known species in their ability to use their tail for rapid and taxing behaviors, including sprinting or climbing.